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Abstract 

We analyze the cross section for vector meson photo-production off a deuteron for the inter- 
CD 

I/"} , mediate range of photon energies starting at a few GeVs above the threshold and higher. We 



reproduce the steps in the derivation of the conventional non-relativistic Glauber expression based 
on an effective diagrammatic method while making corrections for Fermi motion and intermediate 
energy kinematic effects. We show that, for intermediate energy vector meson production, the 



usual Glauber factorization breaks down and we derive corrections to the usual Glauber method 



to linear order in longitudinal nucleon momentum. The purpose of our analysis is to establish 

5-H 



methods for probing interesting physics in the production mechanism for (^-mesons and heavier 
vector mesons. We demonstrate how neglecting the breakdown of Glauber factorization can lead 
to errors in measurements of basic cross sections extracted from nuclear data. 

PACS numbers: 11.80.Fv, ll.80.La 

Keywords: Generalized Eikonal Approximation, Vector Meson Production 



*E-mail rogers@phys.psu.edu 



1 



I. INTRODUCTION 



Coherent vector meson production from nuclei has proven to be a useful tool for studying 
the structure of vector mesons. In the very high energy, small angle scattering regime, well 
above the threshold for vector meson production, the large volume of available experimental 
data involving proton targets consistently supports the validity of the vector meson dom- 
inance (VMD) model for small photon virtualities This, combined with the onset 
of the eikonal regime in the diffractive region has lead to the development of a successful 
theoretical framework for the description of vector meson photo-production off nuclei based 
on the combined VMD model and Glauber theory of hadron-nuclei rescattering . The 
simplicity of the VMD-Glauber framework arises from the fact that at high energies the 
basic 7iV — > VN and VN —>■ NN amplitudes vary slowly with the total energy of the jN 
system relative to the range of important energies in the deuteron wavefunction. This ob- 
servation leads to the factorizability of the basic jN — > VN amplitude from the momentum 
space integral, and yields the conventional Glauber multiple scattering series consisting of 
non-relativistic form factors and elementary scattering amplitudes. 

The VMD-Glauber theory has lead, in particular, to the demonstration that the coher- 
ent photo-production of vector mesons off the deuteron at large —t is defined mainly by 
the rescattering contribution ^j]. Since the VN —>■ VN amplitude appears in the double 
scattering term,one may use nuclear photo-production reactions to study the properties of 
vector mesons p|. By choosing different t, one can control the relative distance at which 
rescattering may occur, which allows one to investigate the space-time evolution of hadronic 
systems produced in electro (photo)- production. 

The above program can be extended to the study of coherent vector electro-production 
at large Q 2 . In this case, coherent vector meson production from the deuteron can be used 
to study color coherence/transparency phenomena in vector meson electro-production at 
high Q 2 . The onset of color transparency will reveal itself through the substantial drop in 
the double scattering contribution with an increase of Q 2 as opposed to the nearly energy 



a. 



independent behavior of the double scattering term for the generalized VMD prediction 

In this paper we consider yet another venue of application for vector meson photo- 
production off nuclear targets by considering photo-production in the intermediate range 
of energies starting a few GeVs above the threshold. These reactions have great poten- 
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tial for probing several effects such as non-diffractive, OZI violating mechanisms for vector 
meson production mesons, in-medium modifications of vector mesons, the importance of 
"non-ideal" u — <f> mixing, and other new mechanisms for vector meson production (see 
Refs. BUS). 

Finally, it would be interesting to learn whether the 0-meson is produced with a small 
enough transverse size that quark degrees of freedom may become relevant, as in the case 
of J/V'-production. Actually, in the case of J/tf), the cross section of the J/if> — N interac- 
tion (Tj/^-N ~ 3mb, [9| estimated based on the A-dependence of J /if) photo-production at 
energies ~ 20 GeV, is much larger than the estimate based on the VDM: < 1 mb. This is 
likely due to the color transparency phenomenon 3- A natural question is whether a trace 
of this effect remains in the case of (^-production. Jefferson Lab has produced data for <f> 
production that is currently being analyzed. 

The interest in intermediate energy reactions makes it necessary to re-evaluate the as- 
sumptions of the traditional Glauber series method, and to develop a new theoretical ap- 
proach. This paper addresses the issues one must face when considering photon energies 
large enough that the eikonal approximation is an appropriate description of hadronic re- 
interactions, but not large enough that it is appropriate to neglect vector meson masses 
in kinematical calculations or any non-trivial s-dependence of the amplitude for photo- 
production of vector mesons from the nucleon. Furthermore, for small photon energies 
(< 3 GeV) the VMD hypothesis becomes suspect as a description of the 'jN — > VN am- 
plitude. Therefore, we will not restrict ourselves to VMD model of 7V — > VN amplitude, 
considering instead the adequately parameterized form of photon-nucleon amplitudes. We 
argue in this paper that there may be a range of photon energies for which the eikonal ap- 
proximation is valid, but where the usual Glauber theory assumptions of factorization and 
ultra-relativistic kinematics break down. 

Although we retain the eikonal approximation, our approach is distinctly different from 
the usual Glauber- VMD approach. In particular, one of the basic assumptions used in 
the Glauber approach is that the basic 'yN — > VN and VN — > VN cross sections are 
slowly varying functions of center of mass energy and that the small Fermi momentum 
of the nucleons can be neglected in the evaluation of the total center of mass energy of 
the 7./V and VN systems. These assumptions result in the usual factorizability already 
discussed above. At intermediate energies, however, the photon energy is comparable to 
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the vector meson mass, and the basic amplitude may gain non-trivial energy dependence 
due to the fact that Regge theory may be inadequate at intermediate photon energies. 
The usual smooth, slow rise in the total 'yN — > VN cross section characteristic of high 
energy diffractive scattering may be absent at intermediate energies. Fermi motion effects 
thereby destroy the factorizability of nuclear scattering into basic amplitudes and form 
factors. Also, the longitudinal momentum transfered (proportional to My/Ej) plays an 
important role as compared to reactions in the diffractive regime and further calls into 
question the factorization assumption. Earlier work (e.g. ^(J) has considered the effect of 
longitudinal momentum transfers, but the breakdown of factorization has not been discussed. 

To summarize, the particular reaction we are interested in in this paper is the coherent 
photoproduction of vector mesons from the deuteron. However, the energy dependence of 
the 7iV — > VN will require that we account for Fermi motion effects which, in turn, will 
require that we account for non-factorization effects. In the derivation of the total 7-D — > VD 
amplitude we will use the Generalized Eikonal Approximation (GEA) with effective Feynman 
diagram rules (see, e.g. This is the approximation, valid at appropriately 



high energies, that allows us to derive the scattering amplitudes starting with corresponding 
effective Feynman diagrams while neglecting multiple scattering from the same nucleon. 
This is very similar to the GEA approach that has been applied to the A(e, e'p)X reactions 
on the nucleus 0|. 

By maintaining the result in terms of momentum space integrals, within the GEA, trans- 
ferred longitudinal momentum and Fermi motion effects may be explicitly taken into account 
consistently. In our derivations, we keep only the corrections to the basic amplitudes that 
are of linear order in longitudinal exchanged momentum or nucleon momentum (neglecting 
order ~k^ 2 /m% corrections, where kN is the bound state nucleon momentum). This allows us 
to relate the D — > NN transition vertex to the nonrelativistic wavefunction of the deuteron. 
Since dynamical, model-dependent corrections related to the N-N interaction are expected 
to be of quadratic or higher order in nucleon momentum, then linear order corrections aris- 
ing from intermediate energy kinematics should be taken into account before any specific 
theory of the basic bound state amplitude that deviates from the nearly flat behavior of 
Regge theory is considered and used in the typical Glauber theory approach. 

As it was explained above, we work in the kinematic regime in which diffractive behavior 
is not yet fully established but the momenta of the produced vector mesons are high enough 



4 



that the eikonal approximation for the hadronic rescatterings is justified. As a result, a 
formalism should be maintained that allows the 7iV — > VN and the VN — > VN amplitudes 
to be independently modeled. Fitting data to our modified form of the Glauber theory by 
using the VN — » VN amplitude as a parameter allows one to infer a value for the VN — > VN 
cross section. We emphasize that the main steps of this paper have been known for several 
decades; the Glauber theory in terms of effective Feynman diagrams was established in 



Ref. 111. The effects of longitudinal momentum transfer in terms of phase shifts have also 

n 

been studied [10||. However, as far as we are aware, there has never been direct numerical 

study of the effect of the breakdown of factorization in Glauber theory as it applies vector 

meson production. (The effects of factorization break-down in proton knock-out have been 

studied in Ref. One result that we find is that the breakdown of factorization persists 

even in the limit that off-shell effects in the bound state nucleon amplitudes are negligible. 

The paper is organized as follows: In Sec. |H] we derive scattering amplitudes for the 

— > VD reaction based on the generalized eikonal approximation. In Sec. IHII we discuss 

the steps needed to take into account linear order corrections in nucleon momentum. In 

Sec. IIVI we perform a sample calculation where we compare our results at intermediate 

energy kinematics with the prediction of conventional VMD-Glauber theory. We identify 

the effects which are responsible for the divergence of the our approach from the standard 

Glauber theory. We demonstrate that effects calculated in this paper, if unaccounted for, can 

yield a misinterpretation of the VN scattering cross section if it is extracted from the data 

using the usual Glauber approximation. In a related issue, we will discuss the recent data on 

fl 

(j) production at SPring-8/LEPS (la) and demonstrate the need to consider kinematic effects 
in the intermediate energy region in Sec. In particular, this data suggests the importance 
of non-vacuum exchanges corresponding to rj and it in the 0-meson production mechanism 
at E 1 of a few GeVs. In Sec. IVII we summarize our results. In Appendix El we describe the 
model of the basic amplitude that we used for our sample calculation, and in Appendix El 
we give an overview of the usual treatment of deuteron spin in Glauber theory. 
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FIG. 1: 3-momentum of a vector meson produced by a photon scattering off a nucleon target as 
a function of — t for a given set of fixed E 1 . The solid lines correspond to meson production, 
whereas the dotted lines correspond to p°-meson production. The incident photon energies in each 
case, going from the bottom curve to the top curve are 1.6,1.8,2.0,2.2,2.4 GeV. Details are discussed 
in the text. 

II. FORMULAE FOR THE AMPLITUDES 
A. Reaction and Kinematics 

We study the coherent photo-production of vector mesons off the deuteron in the reaction: 

7 + D -> V + D', (1) 

where P = (Md, 0) and P' = (Ed, P') define the initial and final four momenta of the 
deuteron. We use natural units (c = h = 1). q = (E 7 , q) and P v = (E v , P v ) define the 4- 
momenta of the initial photon and the final state meson respectively. The three-momentum 
transfered is defined as 1 = q — P V - 

In our calculations we concentrate on intermediate energy kinematics in which, although 
the photon energies are not high enough for the diffractive regime to be established for 
the photoproduction amplitude, the produced vector meson is sufficiently energetic that the 
eikonal approximation can be applied to the calculation of final state hadronic rescatterings. 
This require further elaboration: The GEA is the "straight line" approximation in that the 
incident particle follows a nearly straight line path through the nucleus. Clearly this must 
occur at high enough energies that higher partial waves than just the s-wave contribute. In 
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order to establish the appropriate kinematical regime for our approach, we have plotted in 
Fig. the lab frame 3-momentum of the final state vector meson as a function of —t for a 
set of incident photon energies for the case of p° and 0-meson production from a nucleon. 
Experience with the application of the Glauber model to the description of proton-nucleus 
scattering Q| as well as A(e, e'p)X reactions indicates that the eikonal approximation 
works roughly for Pn/^n ^ 1 {pn is the proton 3-momentum) and it works extremely well 
for En > 2m n- Since ~ rriN, we expect to see the onset of the applicability of the GEA 
for a similar range of momenta for the case of 0-meson production. By analogy with the 
proton case, we continue to use the criterion that E 1 > 2My, and we find that the value of 
vector meson 3-momentum above which the GEA may certainly be applied is Py > 1.8 GeV 
for 0-meson production. The values of 3-momentum, 1 GeV and 1.8 GeV have been indicated 
by horizontal dashed lines in Fig.^ These dashed lines in Fig.^may be viewed as separating 
kinematic configurations where our approach may be applied to 0-meson production from 
kinematic regions where both the approach of this paper and the standard Glauber approach 
should be abandoned entirely with regards to 0-meson production. Below Py ~ 1 GeV, both 
the approach of this paper and the usual Glauber approach should be abandoned. Between 
1 GeV and 1.8 GeV, the GEA may become a rough approximation, but above 1.8 GeV, 
the approach that we take in this paper by using the GEA is a very good approximation. 
For the p°-meson, the eikonal regime begins at smaller values of momentum than for the (f) 
meson due to its smaller mass, so to avoid confusion we do not include the corresponding 
range of applicability of the eikonal approach to p°-meson production in Fig. ^ 

Our main interest in this paper is the production of the 0-meson at around 3 GeV, so the 
application of the GEA is quite safe. We will find that another problem arises at t ~ t m i n , 
and this will be discussed in Sec. IIV1 but the above argument remains applicable as long as 
—t is more than a few tens of MeVs larger than —t m i n . We further assume that the non- 
relativistic model of the N — N interaction can be represented by a D — > NN vertex. In the 
case of the deuteron, there are only two relevant diagrams: the single scattering diagram 
(Born term) of Fig. and the double scattering diagram of Fig. El 
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FIG. 2: The impulse diagram for photo-production. The cross on the spectator nucleon line 
indicates that the spectator nucleon will be taken on shell in the non-relativistic approximation. 
(For all Feynman graphs we use Jaxodraw |ig|.) 
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FIG. 3: The double scattering diagram for photo-production. The crosses on the spectator nucleon 
lines indicate that we will take poles corresponding to these nucleons going on-shell. (See sec. MCjl 

B. The Born Amplitude 

We start with the calculation of the amplitude corresponding to the Born term of Fig. |21 
F^ im ,(s,t) will denote the — > VD scattering amplitude for the Born term in which 
only one of the nucleons takes part in the interaction, whereas F(s,t) will denote the basic 
7 A — > VN scattering amplitude. A hat on a variable indicates that it is associated with the 
7 A — > VN subprocess rather than the process of Eq. The superscript, 0, is meant to 
distinguish the Born term from the double scattering term. The initial and final polarizations 
of the deuteron are denoted by m and m' respectively. Because we consider only intermediate 
energies, the F(s, t) amplitude is not necessarily diffractive and we do not assume the validity 
of the VMD hypothesis. We neglect the spin-flip component of the basic amplitude (i.e. F 
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and F are approximately diagonal is nucleon spin.) The D — > NN vertex is denoted by 
T m . All variables correspond to the labels in the Feynman diagram of Fig. |2]for the single 
scattering (Born) term. The free nucleon mass is denoted by m^. 

By applying effective Feynman rules to the graph in Fig. |2] we obtain the covariant 
scattering amplitude, 



d 4 k Tl,(P -k + l)F(s, t)T m (P - k) 



(2) 

z(2tt) 4 [(P - k + I) 2 -m 2 N + it] [(P - k) 2 -m 2 N + it] [k 2 - m 2 N + it] V ' 
+ (p <-> n). 

(p n) refers to the term in which the neutron and proton are inverted. In the remainder 
of this text, Mandelstam variables that appear within an integral are understood to be 
functions of internal nucleon 4-momentum and the incident photon 4-momentum. 

We proceed with the derivation by estimating the loop integral in Eq.(J2J) up to terms of 
order This approximation allows us to evaluate the integral in Eq.© by keeping only 

m N 

the pole contribution which yields a positive energy for the spectator nucleon. We find, 

d 3 k Tl, (P-k + l)F(s, t)T m (P - k) 



(2tt) 3 2k {(P-k + I) 2 -m% + ie] [(P - k) 2 - m% + te] ( 3 ) 

+(n <-> p). 



(Note that ko = tjin up to correction terms of order h 2 /m 2 N .) We now make use of the 
correspondence between the non-relativistic wave function and the vertex function, 

# Or ,) = -r m (p - k) 

m[ rel) ~ 2^h(2^D(p-ky [ ) 



the form of which is established by the Lippman-Schwinger Equation [20] and by demanding 
that the non-relativistic wavefunction be normalized to unity. Here, —D(P — k) is the 
propagator denominator of the struck nucleon. We write k re ; to indicate that the argument 
of the wave function is the relative 3-momentum of the two nucleons. Using Eq. © with 
Eq. Q and using lab frame kinematics yields, 

< m ,(P 7 ,/) = 2 J d 3 k¥ m ,(k-l/2)F(E~ n k,l)y m (V + (n ^ p). (5) 

We stress, at this point, that Eq. (jSJ) does not yet coincide with the conventional VMD- 
Glauber theory because we have abandoned the usual assumptions that allow us to ignore 



the k-dependence in the basic amplitude, which would normally allow us to factor the basic 
amplitude out of the integral and leave us with the product of the basic amplitude with 
the non-relativistic form factor. For heavier vector mesons (like the 0-meson), the vector 
meson mass may not be negligible, and the s-dependence of the basic amplitude becomes 
non-trivial at intermediate photon energies. 



C. The Double Scattering Amplitude 

Having obtained the Born term in Eq. (j3J), we move on to calculate the double scattering 
term of Fig. El Applying the effective Feynman diagrammatic rules, we obtain 

d* Ps dy s rl,(p + i-ti)F(s,t)F(s,t)r m (P-p s ) 



«(2tt) 4 i(27r) 4 [p2 -m 2 N + ie] [pf - m 2 N + ie] [(P - p s ) 2 - m% + ie] 

+ (p <-> n). 



(6) 



[{P + l- P' s ) 2 -m 2 N + ie] [(q-l + p' s - p s ) 2 — M v + ie] 

Figure 01 and Eq. (JBJ) express the following sequence of events: The incident photon scatters 
from a nucleon with center of mass energy, a/5, producing an intermediate state with invari- 
ant mass, My- The intermediate state propagates through the deuteron before scattering 
from the other nucleon with center of mass energy s. (Bars over variables will indicate that 
they correspond to the secondary scattering.) We neglect fluctuations of the intermediate 
state for the present purposes. Now let us integrate over, p Sj0 and p' s . The integration over 
p S fl is similar to the integration over /c for the Born term of Eq.(j2J). For the p' s integration 
one can choose one of the positive energy poles at Mp + l — ^/m 2 N + (1 — p'J 2 + ie and 
a/ m 2 N + Pg 2 — ie at the upper and lower complex semiplane of p' Q . Note that within the ap- 
proximation in which p'f- ' /m 2 N , I 2 /m 2 N terms are consistantly neglected, the integration over 
either pole will yield the same result. We choose the \/rn 2 N + p^ 2 — ie pole (the poles chosen 
for integration are identified by the crosses shown in Fig. EJ) because this choice reproduces 
the usual Glauber formula in the most direct way. Applying the definition in Eq. we 
recover the formula quoted in ^, 

^m,m' 0^7" 

(2vr) 3 y^p 1 ^ [(q~l + p' s ~ Ps) 2 ~ M 2 + ie] ' 
10 



The (p — > n) term is implicit in these equations. Finally, we put this equation into a 
form that makes the next section slightly more manageable by transforming the variables of 
integration from p' s and p s to p = (p s +p' s )/2 an d k = p' s —p s , and we make the redefinitions, 
k —* k + 1/2, and p — > p + I. The result of these changes is: 

(iVi 3 k^(P + \)F{s,t)F{s, t)%Jp - |) 
(2tt) 3 mjv [(g + fc - |) 2 — My + z'e] U 

+ (p <->• n). 

In Eq. (jHJ), we have given the amplitude a superscript, 1, to distinguish it from the Born 
term. 

We will summarize this section by cleaning up our notation and by writing out the correct 
expressions for the kinematic variables in terms of the integration variables, taking into 
account the variable transformations that were needed to get Eqs. (0) and (jHJ). We explicitly 
expand each expression to linear order in nucleon momentum in the lab frame. Furthermore, 
we assume that nucleon 3-momentum and the exchanged 3-momentum are both small and 
of the same order of magnitude relative to all masses involved. Subscripts a denote Born 
amplitude quantities while subscripts b denote double scattering quantities. The variables 
in each expression are established in the particular diagram under consideration. First, we 
have, 

s a ={{P-k) + qf 

(9) 

= m 2 N + 2E 1 m N + 2E^k z + 0(k 2 ). 

Recalling the variable transformations we made in the double scattering term and noting 
that p, k and 1 are all of the same order of magnitude, we have, 

s b = (q + P -Ps) 2 

( k\ ( 10 ) 
= m 2 N + 2E 1 m N + 2£ 7 ( p z - -± J + 0(p 2 ). 

Note that there is only dependence upon k z and that k contributions come into play only 
at higher order in nucleon momentum. For the rescattering amplitude, we get, 

\2 



Sb = (k v +p 

\- mlr 4- ?.E^m.AT — 9.E., U. — ^ 
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/ k \ ( n ) 
M^ + m 2 N + 2E v m N - 2E 7 [p z - -jf + 0(p 2 ). 



This last value is obtained after the pole in k z is taken, giving the intermediate state an 
invariant mass of k v = M v . The values of t to be used in each of these cases is, 

L = t (12) 

4=g-t) > =g) , + ^-g-* J .) , + OW (!3) 

f '=G +fc ) 2= © 2 -^-ft + ^) 2+ ° (k2) - (14) 

In the usual VMD-Glauber theory expression for the double scattering term, one keeps 
only the perpendicular components of t and t. The terms proportional to k z are small and 
since they come with opposite sign, they tend to cancel if the t-dependence of the basic 
amplitude is nearly exponential. The terms with I 2 , are proportional to t 2 /M D . Thus, we 
continue to neglect both of the first two terms in Eqs. ()13j) and ([14)1 . Finally, we stress that 
(P — k) 2 = (Mo — tun) 2 + 0(k 2 ) so that the struck nucleon may be treated kinematically 
as being on shell up to terms quadratic in the nucleon momentum. 

By using the kinematic expressions of Eqs. l9l through ITU in Eqs. EJandlHl we may ensure 
that the factors multiplying the deuteron wavefunction in each of the integrals is correct to 
linear order in nucleon 3- momentum (or exchanged 3-momentum). 



III. NUMERICAL ESTIMATES AND THE RELATIONSHIP WITH VMD- 
GLAUBER THEORY 



A. Differential Cross Section 



Now that we have calculated the Born and Double scattering amplitudes, let us set up 
notation that allows us to express the total differential cross section in terms of the basic 
amplitudes for 'yN and VN scattering. We do not discuss any physics in this section, 
but simply formulate our notation to allow for convenient comparisons between the present 
approach and the standard Glauber- VMD approach. 

For any exclusive two body reaction involving incoming particles of mass mi and m2 and 
center of mass energy squared, s, the differential cross section may be represented as follows: 

da m ' m ' 1 

— ; — = ^7 -\F mm ,(s,t)\ 2 , (15) 

dt 167T$(S, mi, 1712) 
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where, 

$(s, mi, 1712) = {{s — m,]) + m 2 — 2sm 2 — 2m 1 m 2 ) . (16) 

In particular, the differential cross section for the reaction in Eq. is, 

da m ' m ' 1 

= Tr a>r n ilC'(^)+C'(^)l 2 ( 17 ) 

at l07T < 5P(S, (J, Tfl^jj 

It follows from Eqs. © and (jS)) that the numerical calculations of Eq. ()17j) will require as 
input the amplitudes for for both the 'jN — > VN and the VN — > VN interactions. 

To proceed, we construct a parameterization of the photo-production differential cross 
section in a form that will provide a smooth transition to the VMD-Glauber regime by 
writing, 

, h -yN^VN ~2 / S \ 2(d(i)-l) 

(*.*) = ^:(f) /WM), (18) 



(it 167T yso, 

for the basic 7iV — > VN interaction. In the high energy photon limit, the function f(t) re- 
duces, by construction, to the usual exponential dependence, e Bt ^ 2 , with the constant B that 
is typically used to parameterize experimental data as in, for example, Ref. [jj. The Regge 
trajectory is at(t) = a't + ao. The factor of (;^) a ^~ 1 is the Regge parameterization obtained 
in the VMD-Glauber regime and g(s, i) is a function which adjusts for other s and t depen- 
dence that may appear in the intermediate energy regime, but such that g(s, 0)(^) a ^ ^ 1 
reduces to 1 in the high energy photon limit. By substituting Eq. (JTSj) into Eq. (fT5j) . we 
obtain, 

F^ VN (s,t) =n (s-m 2 N ) (j) f(t)g(sM + v)- ( 19 ) 

The overall normalization is labeled no and is not necessarily related to a total cross section. 
The variable, f), is a possible real contribution to the amplitude. Because P v > 1 GeV for 



the kinematic regime under consideration (see Sec. Ill A|) . the parameterization we use for 
the VN — > VN simply takes a nearly diffractive form, 



F VN ^ VN (s,t) =a VN m + fj)y/$(s, m N , M v )f(s, i). (20) 

The function, f(s,t) reduces by construction to a Regge parameterization, (^) Q (*) _a (°)e' B */ 2 
in the VMD regime. By applying the optical theorem to Eq. (|20p. we see that <tvn(s) is, 
indeed, the total VN cross section. The variable, fj is a possible real part of the amplitude. 
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Our peculiar choice of notation is made so that we may smoothly recover the usual Regge 
parameterizations when we consider the VMD-Glauber approximation. Indeed, applying the 
VMD hypothesis in the appropriate kinematical regime allows us to assume that F(s, t) oc 
F(s,t). Thus, applying the optical theorem would allow one to deduce the VN — > VN 
amplitude. With the standard high energy approximations, we have, 

(21) 

F(s, t) sa VN {i + ffi^e** 

Here, we have put so = 1 GeV for convenience as is often done in parameterizations. In 
this way, we show how our parameterizations reduce smoothly to the expressions obtained 
within Regge theory and the VMD hypothesis. 

One may fit all of the functions that define the expression for F(s, t) directly to data for 
dq/y The function, g(s,i) has been introduced to account for peaks in the energy de- 

pendence or other irregular energy dependence at intermediate energies. Without the VMD 
hypothesis, we can assume no relationship between F(s,i) and F(s,t). At intermediate en- 
ergies, therefore, F(s, t) must be obtained from a theoretical model or by other experimental 
means. Conversely, one can use data for the reaction in Eq. (JTJ) to extract F(s, i). 



B. Corrections to Factorizability and an Effective Form Factor 

We now define an effective form factor, 

stf (*» I) - 1 r5 t^ M fe) ^ w ( k - 8 « k » • (22) 

for the Born term, and an effective basic amplitude, 

F^E^t) = 2E,m N n (i + fj) [^ K J f(t). (23) 

If we substitute Eq. ()19|) into Eq. (j3J), then the Born amplitude for production from the 
deuteron is, 

K,m'( E ^ = 2p ef f (E„ l)S^ a (E 7 , 1/2) + {n~p). (24) 

The definition in Eq. takes the form of a general diffractive parameterization obtained 
when one makes the VMD hypothesis. However, Eq. (|24|) is exactly correct without any 
approximations. We have recovered the usual structure of the Born expression - the product 

14 



of a diffractive basic amplitude with a form factor. The new feature in Eq. (|24jl is that our 
effective form factor depends on the energy of the photon. The definitions that we made in 
Eqs. (J22)) and (|23|) ensure that the effective form factor and the effective diffractive amplitude 
reduce to the usual non-relativistic form factor and the true diffractive basic amplitude in 
the limit that E 7 » M v : 

S™f (£ 7 , 1/2) S m > m ' (1/2) , 

F° ff (E„ I) F VN - VN (s, t). 

By following the usual methods of VMD-Glauber theory, one will extract the effective am- 
plitude from the jD —>■ VD cross section rather than the true amplitude. If, in the region 
of very small —t where the Born cross section dominates, the amplitude for the 'yN — > VN 
scattering is inferred from data using the usual VMD-Glauber theory, then Eq. |^] can be 
used to obtain a corrected amplitude that accounts for non-factorizability. 



C. Corrections to Factorizability in Double Scattering 



The double scattering term is more complicated due to the fact that, in Eq. (JHJ), the 
energy dependence cannot easily be factorized out of the integrand. We may rewrite Eq. (jHJ) 
using Eqs. (JHJ) and (J2H as, 



F^ m ,(E^t) = - / dk z / d 2 k 



d 3 p f(t b )f(h,t b )*Up + |)jUp - |) 

Ml + it] 



x 



27r) 3 m N [(q + k 

a (t*)-i 



1)2 



(26) 



y/$(s, m N , 0)$(s, m N , M v ) 

xg(s, i)n a VN (s) (i + fj)(i + fj). 

The nonfactorizability of Eq. (126(1 near threshold comes from the fact that the basic ampli- 
tudes and the factors in braces have non-trivial dependence upon the integration variables. 
We determine that there is no simple reformulation of the integral in Eq. (|26J) which con- 
sistently accounts for corrections linear in momentum. Therefore, we conclude that a direct 
numerical evaluation is necessary. Note that, though we have set up the integral for a spe- 
cific parameterization, the analysis applies to any smooth, slowly varying energy dependent 
basic amplitude. The k z integral is determined by expanding the denominator in Eq. (|26)). 



15 



(q + k - ^ — My + ie& 2E 1 -k z + l f - ^ + (k - |) + ie 

= 2£ 7 [—k z - A + ie] . (27) 

The second line fixes the definition of A. Notice that by ignoring the term, (k — |) 2 /2i? 7 , 
we have ignored the possibility of contributions from intermediate mesons which are far 
off shell and which correspond to nucleon 3-momenta that are strongly suppressed by the 
deuteron wavefunction. Furthermore, note that the pole value of k z in this approximation 
only depends on the external variables and is independent of the transverse motion of the 
nucleons. The resulting double scattering amplitude is then, 

d 3 P r dk z f(h)f(s b ,mUp+^ m (p-^) 



_! , f , 2l f d 6 p f°° dk z 



x a(t)-l 
S 1 



X | — 



r) 2E 1 m N [k z + A — ie] 

y/$(s, m N , 0)$(s, m N , M v )g{s, t)n Q a VN (s b ) (i + fj) (i + fj) ( 28 ) 

dk z / m , m >(k ± ,fc z ,p,s) 
(27r) 2E 1 m N [k z + A-ie]' 

We have gathered all factors apart from the energy denominators in the integrand into a 
function, / m , m '(kj_, k z , p, s). Assuming identical protons and neutrons, we get an identical 
term for the case where the roles of the neutron and proton are inverted. A convenient way 
to reorganize this formula so that it more closely resembles the non-relativistic quantum 
mechanical theory is to write the integrand in terms of its Fourier components in the following 
mixed representation: 

1 f°° 

Im,m>(k±, k g , p, s) = — = / cb/ m , m /(k ± , z, p, £ 7 )e" ifc ^. (29) 

V 27T J-oo 

The vector meson propagator may be rewritten using the identity, 

1 f°° 

= / dzQ(-z)e^ p - ie >. (30) 

p-ie y„ 00 

Summing the two terms for the neutron and the proton and using the fact that Q(z) + 
Q(-z) = 1 yields, 

<i 3 p<i 2 kj 



i f d pd k ± . 

F m , ml {E^) =zj ^_^/ m>m ,(k J .,-A,p, J B 7 ) 

1 f d 3 pd 2 k ± f°° , - „ , . , . . 

~ V2l J 2 El m N (2nf L dzI ™>™'^ *> *>> S ) ^(-A,)0(,). 



(31) 
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In the VMD-Glauber approximation, A — > and I m , m '(kj_, 0, p, s) is the usual energy inde- 
pendent density matrix. Hence, the first term in Eq. ()31j) reduces to the traditional Glauber 
expression for double scattering and the second term vanishes in the limit where the usual 
VMD-Glauber assumptions are applicable. The second term is a correction, discussed in 
Ref. P| which arises from the non-zero phase shift in the vector meson wave function induced 
by longitudinal momentum transfer. In the phase shift term, the factor of sin(— Az) is itself 
a correction of order k z , so we neglect Fermi motion and energy dependent corrections to 
I m , m '(kj_, 0, p, s) in the phase shift term. 

For a real photon, the double scattering term picks out the relative longitudinal nucleon 
momentum, 

A = 

Furthermore, I- = — {M v — t)/2E 7 , so 

A 

and we see that A is indeed negligible at large center of mass energies and small t. Correc- 
tions to the double scattering term, at linear order in momentum, arise from performing the 
integral in the first term of Eq. (}3"T]) numerically, and by retaining the phase shift term. We 
end this section by noting that the breakdown in factorization comes simultaneously from 
the fact that longitudinal momentum transfer is non-negligible, and the fact that the lon- 
gitudinal momentum of the bound nucleons is non-negligible; the contribution to the basic 
amplitudes from the longitudinal component of the bound nucleon momentum at linear or- 
der would vanish by symmetry in all of the integrals if the longitudinal momentum transfer 
were neglected in the wavef unctions. 

IV. SAMPLE CALCULATIONS 
A. Cross Section Calculation 

It is usually the case that one calculates the charge and quadrapole form factors in the 
coordinate space formulation of the form factor. This method reduces the formulae to an 
extremely simple form and allows one to deal simply and directly with polarizations. For 
the purpose of modifying the basic amplitude, however, so that it has nucleon momentum 



~2 



2EL 



t + My 

4£L 



(32) 



(33) 
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dependence, we must maintain the momentum space formulation that results from a direct 
evaluation of the effective Feynman diagrams in Figs. (I2I3D . Carrying this procedure out 
was the topic of the previous two sections. The calculation is straightforward, but becomes 
numerically cumbersome, and the longitudinal momentum exchanged leads to a breakdown 
of the orthogonality relations for spherical harmonics that usually lead to a very simple 
coordinate space formulation. However, dealing with the deuteron polarizations can still 
be simplified if one chooses the axis of quantization along the direction of momentum ex- 



change 3|,|2l|. An overview of the non-relativistic deuteron wave function with polarizations 
are described in Appendix iBl 

In this section we provide some sample calculations by using simple models of the basic 
amplitudes. To this end, we restore the assumption of VMD and we use very simple param- 
eterizations of the s and t dependence in the basic amplitudes. The purpose for doing this 
is mainly to provide estimates of the sensitivity to non-factorizability rather than because 
VMD is thought to be appropriate at intermediate energies. We have extracted estimates 
of the parameters for production of the p° and (j) vector mesons from the basic nucleon in- 
teraction cross section data appearing in Ref. Q], and we have made rough estimates of the 
parameterization of the s and t-dependence of these amplitudes (see Sec.[X]for a description 
of our parameterizations). This provides us with a reasonable model to work with, though 
we stress that refinements are ultimately needed. For all of our calculations we use the 
non-relativistic wavefunction obtained from the Paris N-N potential J3] . 

We are mainly interested in the ^-production cross section which is dominated by natural 
parity exchange, even at energies close to threshold, due to the OZI rule. However, to 
demonstrate the consistency of our approach with traditional methods, we consider first 
the case of the photo-production of p°-mesons which has been well understood for some 
time. The basic amplitude for p°-production is dominated by soft Pomeron exchange at 
large energies, so that it is constant at high energies, but undergoes a relatively steep rise at 
energies near threshold due to meson exchanges. The parameterization we use is shown in 
appendix |XJ We use a typical exponential slope factor of 7.0 GeV -2 for the t-dependence. 
Fig. (@J) shows the cross section for p°-production at the high energy of E~ = 12.0 GeV. For 



^ we show data ta.ea at !,0 GeV fr0m Re, fl M . The co mparis o„ with data is 
reasonable, as it is with the traditional Glauber approach. 

Now we consider the more interesting case of 0-meson photo-production. At high energy, 
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Total Cross Section (p°) 




10" 

0.2 0.4 0.6 

-t (GeV 2 ) 

FIG. 4: The unpolarized differential cross section for coherent /3°-meson production compared with 
the total cross section for different polarizations. The calculation is done with the large photon 





energy E~ = 12 GeV, and the data for E~ = 12 GeV is taken from Ref. 



we use the Regge dependence, a(t) = .27 1 + 1.14 given in Refjjj. The parameterization that 
we used is described further in appendix|Xl As noted in Ref. J|, the energy dependence of the 
0-meson photo-production cross section is very weak, but the current state of experimental 
data is still ambiguous as to how much this energy dependence continues at lower energies. 
However, the large negative ratio of the real to imaginary part of the amplitude (rj = 
—.48) 23[ suggests that some mechanism other than soft Pomeron exchange is significant. 
This value of the ratio of the real to imaginary part of the 0-meson cross section has large 
error bars and was calculated neglecting longitudinal momentum transfer. However, it is the 
only measurement we know of at the moment so we use it for the purpose of demonstration. 
At lower photon energies than what we consider here, the energy dependence of the basic 
cross section may become highly non-trivial as is suggested by data in Ref. . The results 
of the calculation done with each combination of initial and final deuteron polarizations are 
shown in the separate panels for a photon energy of E 1 = 30.0 GeV in Fig. El and for a 
photon energy of E 1 = 3.0 GeV in Fig. |H1 The result is summarized in Figs. [3 and El which 
show the differential cross section for different polarizations along with the unpolarized cross 
section for photon energies of 30.0 GeV and 3.0 GeV respectively. 
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Each of the curves in Figs. El and El separately represents the contribution to the total 
cross section from a term in the squared amplitude when we apply Eq. El The Born and 
double scattering terms are obtained from the square of Eq. and square of the first term 
of Eq. |inj respectively. The phase shift term arises from the square of the second term in 
Eq. |^ We call it the phase shift term because, in the language of non-relativistic quantum 
mechanical wavefunctions, it arises due to a phase difference between the incoming photon 
and the produced vector meson. The interference term arises from the interference between 
Eq. |^ and Eq. |^ Note that the interference term is negative, but it is plotted on the 
positive axis for demonstration purposes. Note also that there is no contribution from the 
Born term for the m = +/-ltom = -/ + l transition, and therefore the total cross section 
for the spin-flip reaction has none of the large dips characteristic of the Born cross section. 

An important feature that can be seen in Figs. El and El is that the double scattering term 
is suppressed in the intermediate energy case relative to the high energy case. We can see this 
most clearly by comparing the upper left panel of Fig. El with the upper left panel of Fig. El 
It is clear that the double scattering contribution is important in the = 30 GeV case at 
moderate values of —t, whereas for the E^ = 3 GeV case the cross section is dominated by 
the Born term all the way up to —t ~ .4 GeV 2 . In the general case of multiple scattering from 
complex nuclei, it is the rescattering contributions which lead to the usual A-dependence (A 
is the number of nucleons) of Glauber theory. The fact that multiple scattering is suppressed 
in double scattering in the deuteron suggests that our method would yield a rather different 
A-dependence from that of usual Glauber theory if it were extended to complex nuclei. 
Extending our approach to complex nuclei will be the subject of future work. 

Another problem begins to emerge at lower photon energies and extremely small —t (at 
t ~ tmin)'- A large fraction of the momentum integrals begins to violate relativistic kinematic 
constraints. It is likely that the basic amplitudes vary extremely rapidly with s and —t in 
these regions of the integral and that expanding in nucleon momentum is not valid (at least to 
linear order). In order to make progress, a precise understanding of the dynamics of off-shell 
amplitudes based on field theory may be necessary. Therefore, our approximation is only 
valid at —t sufficiently large that the integrand does not contain significant contributions 
from kinematically forbidden nucleon configurations. We have tested the effect of this region 
in our calculations, and in performing our calculation, we find that there is virtually no 
contribution from kinematically forbidden regions for any situation that we consider as long 
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as — t + t m i n is greater than a few tens of Me Vs. We note that, even at relatively low photon 
energies, the data is consistent with a smooth exponential — t-dependence (see Ref. ) 
as long as —t is not exactly — t m j n . Note that this theoretical problem of considering t at 
exactly —t m i n exists at high energies as well, but that at high energies —t m i n is generally 
too small for it to show up in plots. So that we may perform our calculations numerically 
at all values of —t greater than —t min , we choose to make the basic amplitude vanish in 
kinematically forbidden configurations (when —t < —t m i n ). This results in a small dip 
just above —t m i n in our plots. The small dip is, therefore, unphysical, and should not be 
regarded as a prediction. We leave it in our plots merely to illustrate a general failure of the 
Glauber theory approach at extremely small —t (see Fig. [H] at —t < .06 GeV 2 ). In summary 
of the above, the small dip at extremely small —t denotes a kinematic region in which no 
known multiple scattering formalism works. Numerically, our calculation is only correct in 
the region of —t above the dip at small —t; that is, when —t + t m i n is greater than a few 
tens of Me Vs. 
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m = +/-1 to m = +/-1 




0.2 0.4 0.6 0.2 0.4 0.6 

-t (GeV 2 ) -t (GeV 2 ) 

FIG. 5: The long-dashed, dotted, dashed, solid, and dot-dashed lines refer to the Born, double, 
interference, total, and phase shift terms respectively for a photon energy of E 1 = 30.0 GeV. The 
interference term is negative but is plotted for illustration on the positive axis. See text for detailed 
discussion. 
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FIG. 6: The long-dashed, dotted, dashed, solid, and dot-dashed lines refer to the Born, double, 
interference, total, and phase shift terms respectively for a photon energy of, E 1 = 3.0 GeV. Note 
the different scale on the axis for the spin-flip contribution. The interference term is negative but 
is plotted for illustration on the positive axis. See text for detailed discussion. 
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(^-meson photo-production 
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FIG. 7: The differential cross section for <^>-meson production for different polarizations for a 
photon energy of E 1 = 30.0 GeV. See text for detailed discussion. 
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(^-meson photo-production 




0.2 0.4 0.6 

-t (GeV 2 ) 

FIG. 8: The differential cross section for 0-meson production for different polarizations For a 
photon energy of E 1 = 3.0 GeV. The deuteron spin flip term is negligible at these energies. See 
text for detailed discussion. 
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Next we plot consider the total unpolarized cross sections as a function of photon energy 
for a set of fixed values of —t. This allows us to compare the factorized and unfactorized 
calculations directly and to determine at approximately what value of energy the transition 
to the VMD-Glauber regime occurs. Recall that it is the motion of the nucleons in the 
deuteron (the Fermi motion) that leads to the non-factorizability of the basic amplitudes. 
Factorization refers the practice of ignoring the dependence of nucleon momentum inside 
the basic amplitudes when integrals over nucleon momentum are performed. The use of 
non-factorized amplitudes is the essential difference between our approach and the usual 
Glauber approach. The ratio of the cross section with the usual factorization assumption 
to the cross section which accounts for non-factorizability effects (Fermi motion) is shown 
in Figs. IH1 The upper panel refers to the case of p° production, whereas the lower panel 
refers to ^-production. The ratio is given for two small values of —t: t = — .04 GeV 2 and 
t = -.14 GeV 2 . 

The upper panel demonstrates that the effect of non-factorizability is small for the case 
of p° for the entire range of intermediate energies. This is in sharp contrast to the case of </> 
production in the lower panel of Fig. El Note that we only plot the case of t = —.04 GeV down 
to E y = 5 GeV for the 0-meson case. This is because, for photon energies lower than 5 GeV, 
t = —.04 GeV becomes too close to t m i n . On the other hand, for the curve corresponding to 
t — — .14 GeV, there is nearly a 30 percent suppression of the factorized cross section relative 
to the unfactorized cross section at the lowest energy, E 1 = 3 GeV, shown in the lower panel 
of Fig. El We emphasize that this result is for a photon energy (3 GeV) that is well into 
the kinematic region where the eikonal approximation may be applied (see Sec. Ill Ajl . and 
that —t = .14 GeV 2 is certainly large enough relative to —t min = .036 GeV 2 that there are 
none of the problems discussed earlier related to nearness to —t m i n . Therefore, our method 
of calculation is ideally suited to the kinematics of the dotted curve in Fig. El where a 
significant effect from the break down of the factorizability assumption is already seen. 

Note from the general behavior in Fig.EJthat the cross section rises when the factorization 
assumption is removed. This effect is mainly due to the suppression of multiple scattering 
when non-factorizability is taken into account. To see this, note that At —t = .14 GeV 2 
the main effect of double scattering in the usual Glauber approach is to produce a large 
negative cross term that has a canceling effect. All terms apart from the Born term and the 
interference term are negligible at this value of t. (See, for example, the upper left panels of 
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FIG. 9: The energy dependence of the ratio between the differential cross section calculated using 
the usual factorization assumption and the differential cross section calculated with factorization 
break-down taken into account. 

Figs. El and EJ) Therefore, if multiple scattering is suppressed, as it is in our approach, then 
the absolute value of the cross term becomes smaller, and the Born contribution is no longer 
suppressed by multiple scattering. Thus, the curve representing our approach in Fig. |U] is 
smaller than what is found in the standard Glauber calculation. 

At high energies we expect the two methods to agree, and they do within the range of 
experimental uncertainties of non-relativistic deuteron form factors. The fact that the two 
methods have slight disagreement at high energies is a reflection of the fact that, even at high 
energies, we have not calculated the form factor with exactly the same approximation as in 
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the usual Glauber approach. In the usual non-relativistic form factor, any dependence on 
longitudinal transferred momentum is ignored. If one takes into account exact kinematics, 
one finds that there are two distinct effects which may cause this assumption to be violated. 
It is easiest to see this by writing out the exact expression for the transferred longitudinal 
momentum: 

_t_ u^t 

From this we see that there are two approximations that are normally made in the Glauber 
approach that allow one to neglect l z . The first is the ultra- relativistic approximation for 
the incident vector meson, E 7 » My, and the second is the non-relativistic approximation 
for the exchanged 4-momentum, —t << Mp. If t is small relative to Mp then there is still 
a significant contribution to l z when My is non-negligible relative to E„ r This is the effect 
that interests us in this paper. It is safe to use the non-relativistic form factor because the 
transferred energy is, 

'» = -2^' (35) 
which is small at small —t. On the other hand, as long as t is not exactly zero, there will be 
a component of l z that does not die out with energy. This effect represents the error induced 
by ignoring relativistic recoil. In the future we plan to generalize the formalism to the case 
of light-cone wave functions so that it may be extended to higher —t. 



V. DIRECTIONS FOR FUTURE WORK 
A. Extraction of Basic Amplitudes 

We emphasize that the work in this paper is a first step in refinements to the usual 
techniques applied to multiple scattering in vector meson production from the deuteron. We 
plan to extend these refinements in the future to include, for example, light-cone kinematics 
in the treatment of the deuteron wavefunction, and spin-flip effects. Obtaining precise 
parameterizations of the s and t dependence is one of several steps needed for refinements 
in the calculation. We note that a peak in the energy dependence has been reported in 
Ref. Q| for photo-production of 0-mesons from a proton target at E 1 = 2 GeV and it is 
this data that we used in our parameterization (see Appendix EJ • We would like to point 
out, however, that the measurements in Ref. [la] are for the differential cross section at 
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t = t m i n . Therefore, since the value of t m i n varies significantly with energy in these near 
threshold measurements, then the reported measurements give the differential cross section 
at very different values of t. We have indicated this in Fig. ITUT A.). In order to infer the 
energy dependence at a fixed value of t, one needs to assume a form for the t-dependence. 
The actual t-dependence at these low energies is not well known, but it is straightforward to 
see that even a simple exponential t-dependence will have an effect on the shape of the over- 
all energy dependence of the cross section. As an example, we have plotted in Fig. [TUT A.) 
the data as it was originally presented in Ref. [r| alongside Fig. HDlf B.) where the data 
have been shifted to a fixed value of t. We have used an exponential slope parameter of 
4 GeV 2 which gives reasonable agreement with the data. In the original form of the plot, 
Fig. (A.), the data is shown at a different value of t at each energy. The highest value of 
—tmin occurs at the lowest energy plotted which is around 1.6 GeV. In Fig. HUT BO we have 
re-plotted the energy dependence, but with the value of t for each data point fixed at —tmin 
for Ej = 1.6 GeV since this is the largest value of —t that is kinematically allowed for every 
point on the plot. Let t min [1.6] represent the value of —tmin at E 1 = 1.6 GeV. Then, 



da 
~dt 



da 
~dt 



e 4.0GeV- 2 (t min [1.6]-i mi „) ^g-j 

i — tmin [1-6] t — t Tn { Tl 

We use this to obtain Fig. El (B.). We see that much of the peak-like behavior is removed. 
Without a fuller understanding of the t-dependence, therefore, it cannot be ruled out that 
the observed peak arises from purely kinematical effects. However, the fact that the cross 
section at fixed t does increase at smaller E 1 is evidence that OZI-violating meson exchange 
effects become important at these energies. 



Recently, preliminary data were reported from SPring-8/LEPS |24| which measured the 
dependence of the 0-meson production cross section at E 1 ~ 2 GeV on the linear polarization 
of the photon. Significant polarization is observed which requires the presence of a non- 
vacuum exchange like n, r\ exchange. Such exchanges lead to spin flip in the nucleon vertex. 
These contributions for small t are strongly suppressed for coherent production off the 
deuteron (pion exchange does not contribute in any case due to the zero isospin of the 
deuteron.) These effects are determined by the deuteron magnetic form factor which is 
much smaller than the electric form factor. Hence, the coherent production of the 0-meson 
may be used as a spin analyzer of the elementary amplitude in the kinematics where double 
scattering is a small correction. This topic has already been discussed in Ref. Q. In the 
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FIG. 10: Plot of recent data from LEPS, taken from Ref. 16]. We indicate the significant variation 
of t m in with photon energy. This may have an effect on the over all energy dependence of the cross 
section. (A.) shows how the data was originally presented: at a different value of t for each energy. 
In (B.) we have shifted all of the data points to the same value of t by assuming a constant slope 
parameter of 4 GeV -2 . Each point in (B.) corresponds to the differential cross section at the fixed 
value of t corresponding to to t m i n for a 1.6 GeV photon. Note the different scales on the axes in 
(B.) 

spirit of the original Glauber approach, we have neglected spin effects in this paper for the 
sake of simplicity. Future work will involve generalizations of our method to the case of spin 
dependent basic amplitudes. However, if one fits a combination of Pomeron trajectory and 
Reggeon trajectory to the recent preliminary SPring-8/LEPS data, and then extrapolates to 
3 GeV, then it appears that less than 20% of the basic 7-D cross section is due to spin-flip, 
whereas the corrections found in this paper due to non-factorizability are as large as 30% at 



3 GeV 



24 1 . 



Before ending our analysis, we mention that, because V — N cross sections are extracted 
from the multiple scattering term, quantities sensitive to the deuteron polarization would be 
ideal for testing whether the V — N cross section is unusually large. In order to emphasize 
this, the cross section for scattering from a polarized deuteron, from m = +/ - 1 to m = 
+/ — 1, is plotted in Fig. [TT] where the result of using a typical value for the total 4>N 
cross section (11 mb) is compared with the case when the <pN cross section is enhanced 
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<|>-meson photo-production 

E 7 = 3 GeV, m=+/-1 to m=+/-1 




0.2 0.4 0.6 

-t (GeV 2 ) 



FIG. 11: The energy dependence of the m = +/ — 1 to m = +/ - 1 differential cross section for 
0-meson photo-production with = 3 GeV. The dashed curve shows the result of increasing the 
typical basic 4>N cross section by a factor of 3. 

by a factor of three. (For clarity we have only plotted the sum of all the terms from the 
squared amplitude rather than each term separately.) The sharp dip that normally appears, 
is due to the sharp dip in the Born cross section. However, the double scattering cross 
section is nearly flat in —t. Therefore, in the summed cross section the double scattering 
term dominates in the region of the dip, and may even cause the dip to vanish entirely if it 
becomes very large. Figure^Jshows that, even with the suppression of the double scattering 
term that results from the non-factorizability that we have been discussing, the dip in the 
cross section is observed to flatten out when the basic <pN cross section is abnormally large. 

B. The Problem of Bound State Amplitudes 

We have treated the struck nucleon as being on-shell which is consistent with the neglect 
of terms quadratic in nucleon momentum. However, immediately at the threshold for par- 
ticle production, the 7iV — > VN amplitude has very unpredictable behavior which may be 
modified significantly when the nucleon is in a bound state. This is especially clear when 
we realize that for a given photon energy —t min is different for a deuteron and an on-shell 
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nucleoli target. We cannot predict the effects of the off-shellness of the bound nucleon with- 
out a complete, relativistic understanding of the basic amplitude. However, we have made 
predictions in the region of kinematics where it is reasonable to assume that the bound state 
amplitude is the same as that of the free nucleon amplitude. If one includes dependence 
upon the nucleon virtuality in the basic amplitude, then one may write the amplitude as 
F(s, i, k 2 N ). For k 2 N = m 2 N , the amplitude reduces to the free nucleon amplitude. As we have 
stated, k 2 N = m 2 N up to corrections of order / m 2 N or higher whereas s has linear order 
corrections in nucleon momentum. Thus, if F(s, i, k 2 N ) is an analytic function of kinematic 
variables, then there will be linear order corrections in nucleon momentum due to s whereas 
the lowest order corrections due to the virtuality of the bound nucleon are only of quadratic 
order in nucleon momentum. In other words, Fermi motion effects may be important even 
when it is appropriate to neglect the off-shellness of the bound state amplitude. Of course, 
all of this depends on the validity of using kjv/mjv as a small expansion parameter which is 
only true if the basic amplitude has relatively weak s-dependence. This is one reason why 
we emphasize that we are considering intermediate energies rather than low energies. One 
may also include the deuteron binding energy in the calculation of the mass of the bound 
nucleons, but the binding energy arises from the full consideration of relativistic binding 
and higher order terms in nucleon momentum, so considering the nucleon binding energy 
is not consistent with the neglect of higher order nucleon momentum terms or the use of a 
non-relativistic potential for the N-N interaction. 

In this subsection we propose a rough a way to test the validity of the on-shell amplitude 
approximation. We do this in the next few paragraphs by directly comparing the amplitude 
when it is evaluated at the value of t min for the deuteron with the case when it is evaluated 
at t m in for a free nucleon with s given by the exact expression for s, 

s = 2£ 7 (m d - ^m 2 N + k 2 + k)j + ^M D - ^m 2 N + k 2 ^j - k 2 . (37) 

The value of the nucleon 3-momentum thus parameterizes the off-shellness of the bound 
nucleons. t' min will denote the lower bound of —t for the free nucleon, whereas t min is the 
lower bound of —t for the deuteron. The struck nucleon inside the deuteron for the unprimed 
case has s given by Eq. ()37j) . We will now consider the case of a free nucleon, with the same 
s as for the bound nucleon, but with the nucleon on-shell (i.e. k 2 = m 2 N ) and with a fixed 
value for k z . So that the free nucleon energy corresponds to the bound nucleon energy, we 
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will continue to use Md — a/ rn? N + k 2 for the energy of the struck nucleon. In short, we are 
comparing t m i n for 7 scattering off a deuteron at rest with t' min for 7 scattering off a free 
nucleon, with energy corresponding to that of the bound nucleon in both cases. 

We expect the rate of variation of the basic amplitude with t to be very large near t m i n . 
If there is a significant contribution to the integral in Eq. (jSJ) from regions near t m i n , then 
t m i n should nearly equal t' min in order to make the on-shell amplitude a valid approximation 
to the bound state amplitude. We can use the difference between these two values of t min 
to estimate the effect of the the off-shellness on the the amplitude. 

In order to test the effect of the off-shellness of the basic amplitude, we may consider two 
extremes. First, the bound state basic amplitude could be evaluated at the physical value 
of t for the photon-deuteron process. That is, we could calculate the amplitude, F(s, t) 
at t where t is the physical value of t for the photon deuteron process. In this case, since 
t min is smaller for the deuteron than t' min is for the nucleon, then we are probably over- 
estimating the cross-section. On the other hand, we could evaluate the basic amplitude at 
F(s, t— {tmin— t'min)) = F(s, t—At) where t' min is the minimum t for the free, on-shell nucleon. 
With this second method for choosing which value of t to use in the basic amplitude, the 
basic amplitude behaves like the free, on-shell nucleon amplitude in the region of t close 
to t m i n . Hence, with this method, we are probably underestimating the value of the basic 
amplitude. In the high energy limit, At vanishes and the two amplitudes are equal, and 
the difference between the two provides an estimate of the off-shell effects. (Note that we 
must specify a value for k in order to make a comparison.) Any amplitude which has a 
relatively slow and smooth variation with t will yield a small difference between F(s, t) and 
F(s,t — At). We made this comparison for 0-meson production with the parameterization 
in appendix [X] and we find only a few percent deviation. We conclude that at a few GeV 
above threshold it is reasonable to continue using the on-shell amplitude of the nucleon. 



VI. SUMMARY AND CONCLUSIONS 



The main conclusion of this paper is that the effect of factorization break-down is signif- 
icant for intermediate photon energies. The Glauber approach is, strictly, only applicable 
for the case of very high photon energies. However, there are current attempts to apply the 
factorization assumption of Glauber theory to the 0-meson production reaction at energies 
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as low as 1.5 GeV in both experimental and theoretical research. Therefore, in order to 
salvage the situation in the energy range of a few GeVs above threshold, we have outlined 
steps one must follow in order to obtain corrections to leading order in the bound nucleon 
momentum and transferred momentum. The main steps are essentially those of the original 
diagrammatic formulation of Glauber theory in terms of momentum space integrals and its 



extension to vector meson production 



10 



11 



we have started with most of the original 



assumptions, but we have removed the assumptions of factorizability, ultra-relativistic kine- 
matics or VMD for the basic amplitudes, and we have numerically evaluated all integrals 
directly without any factorization approximations. By using a simple model for the basic 
amplitude (we restore VMD for the simple model) based on a fit to old and recent data, 
we have shown that, away from t = t m i n , ignoring Fermi motion (and the resulting break- 
down of factorizability) can lead to a significant error in basic cross sections extracted from 
— > <pD cross section data (see Fig. (HJ). This effect will certainly need to be taken into 
account in future searches for new production mechanisms at intermediate energies. The 
breakdown in factorizability arises as a consequence of both the non-negligible longitudinal 
momentum exchanged, and the non-negligible Fermi motion. An important point is that a 
key source of the departure from factorizability is the inadequacy of assuming the nearly 
flat s-dependence predicted by Regge theory in the basic amplitudes. Therefore, models of 
the basic jN — > VN amplitude or the VN —>■ VN amplitude which depart significantly 
from nearly flat s-dependence must include at least the linear order nucleon momentum 
corrections of this paper if they are to be used in calculations with a deuteron target. This 
correction arises purely from the fact that the bound nucleons have non- vanishing momen- 
tum and it must be included regardless of the details of a particular model of the basic 
amplitudes. For the case of 0-meson production, we find that our approach is reasonable 
when we use our particular simple model of the basic amplitude and as long as the photon 
energy is around 3 GeV or higher and t is not too close to t min . 

However, we stress that in a model of the basic 'yN — * VN amplitude that predicts much 
wilder energy dependence at intermediate energies than what we have assumed, the linear 
order corrections will not be sufficient, and a complete and precise understanding of the 
N — N interaction and the bound state nucleon amplitudes are necessary in order to make a 
correct calculation. For the 0- meson photo-production cross section (M^ « 1.02 GeV), the 
amplitude may vary wildly with energy at E 1 = 2 GeV or lower because of the very close 
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proximity to threshold. For this reason, and because the eikonal approximation begins to 
break down, basic cross sections for photo-production from the nucleon extracted from data 
for photo-production from the deuteron are suspect for photon energies less than or equal 
to 2 GeV for the production reaction. 

In our sample calculation, we observe that the contribution from double scattering be- 
comes numerically suppressed relative to the Born approximation as the incident photon 
energy decreases. However, the multiple scattering terms are what lead to the characteristic 
A-dependence of the Glauber theory for complex nuclei, a tot ~ A 2 / 3 . This suggests that an 
extension of our methods to complex nuclei will yield a rather different A-dependence for 
the cross section at intermediate energies from what is predicted at high energies. Hence, 
there will need to be a revision in efforts to extract basic cross sections from nuclear data 
using extrapolations in A. The extension to complex nuclei, however, requires much more 
work. We note, however, that data given in Ref. were interpreted as implying a very 
high <f)N total cross section on the basis of a very traditional Glauber approach at energies 
of only a few GeVs. Therefore, our next step will be to determine how the non-factorization 
effects discussed in this paper affect a general, incoherent Glauber series. Furthermore, since 
it is apparent that spin effects will be important, then a generalization with spin-dependent 
amplitudes will be needed. 

We have purposefully over-simplified our analysis here for the purposes of demonstration. 
In particular, we have applied the VMD hypothesis at energies where it is suspect and we 
have neglected fluctuations and u — <fi mixing in the intermediate vector meson in the double 
scattering term. Further analysis will need to include these effects. In order to make further 
numerical progress, we will need firmer parameterizations of the basic cross sections for 
vector meson production from nucleons. For theoretical work, it would be useful for the 
purposes of comparison to have a widely agreed upon set of parameterizations. We also 
need to consider the calculation in light-cone coordinates and the effects of spin-flip. We 
will pursue these issues in future work. 

Finally, we will need a complete understanding of the off-shell amplitudes if we are to 
take into account the higher order momentum corrections that will be necessary just at the 
threshold, though we have argued that for smoothly varying basic cross sections, the effect 
of off-shellness in the amplitudes is small relative to the effect of linear order corrections in 
nucleon momentum. 
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APPENDIX A: PARAMETERIZATIONS 

Here we describe the fits of the basic cross sections that we used for our sample cal- 
culations. The object here is not necessarily to produce very accurate parameterizations, 
but rather to devise parameterizations that demonstrate the effects of Glauber factorization 
while being consistent with recent and established experimental results. 

First we consider the — > p°N differential cross section. For this we use a simple 
exponential t-dependence with a typical exponential slope of B = 7.0 GeV~ 2 and an over 
all normalization of 105 pb/GeV 2 . (See, e.g. Ref. It is known that at low energies 

the normalization undergoes a steep rise. We take this into account in our calculation by 
including a factor of (1 + in the overall normalization and then doing a least squares fit 
to obtain the parameter, a. We find that a rs 32.7. The cross section is thus, 



The result is shown in Fig. As is seen in the main part of the text, the variation is 
too weak to introduce a very large effect on the final 7Z) — > p°D cross section from Fermi 
motion. 

The case of the 0-meson is more complicated due to the irregular behavior near threshold. 
The main point is to interpolate smoothly between recent low energy data and the standard 
higher energy parameterization. The normalization of the low energy data, taken from 
recent experimental work in Ref. |16], is obtained from an effective Pomeron and pseudo- 

n n 

scalar exchange model |26J as it was presented in Ref. [161 ] . We continue to use this so 
that our model will be consistent with current work. The high energy parameterization was 
obtained in Ref. [l] by fitting a diffraction-like cross section to a large set of experimental 
data. We want to interpolate quickly but smoothly between the low energy and high energy 
data. There is an exponential factor, e Bt , associated with both the high and the low energy 
behavior, but the slope, B, is around 3.4 GeV -2 for the low energy behavior (E 1 < 4.0 GeV) 




(Al) 
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while it is around 4.8 GeV 2 for the high energy behavior. Thus, for the exponential slope 
we use, 

£?(£ 7 ) = (4.8 - (4.8 - 3 .4)e-°' 001Gey " 4 ^) GeV- 2 . (A2) 

Next, for the low energy region, there is no Regge slope. That is, a' = in the factor, s a ' 1 . 
But, in the high energy region, a 1 = .27 GeV -2 . Thus, we use, 

a'(E,) = .27(1 - e-ooioev-*^) Ge y-2 (A3) 

Now we consider the behavior of ^- t -_ for photo-production of the 0-meson from a proton 
target. The high energy parameterization in Ref. |l| is 



da 
~dt 



t=o 

and, as in Ref. LU, over- all units will be understood to be G ^ 2 . In order to match to the 



1.34s- 28 , (A4) 

fib 



data of Ref. we want a peak to appear at around E 1 = 2.0 GeV. Therefore, we adjust 



the parameterization to, 

^ =1.34s- 28 (l + ae- b ^- c > 2 ). (A5) 

We use Eq. (IA5|) to fit to the low energy data of Ref. 16] while assuring that the high 
energy parameterization of Eq. (jA4|) is reproduced at high energies. We find: a = .71, 
b = 16.5 GeV" 2 , and c = 2 GeV -2 . Finally, we note that the low energy data is actually 
given for t = t m i n rather than t = 0. Therefore, we must be sure to include a factor of 
e B(E~,)t min j n final result for ^|£_ - The result of our parameterization for ^-| t -_ for 
the 0-meson is shown in Fig. We point out that in the intermediate energy range at 
around E 1 = 3 GeV, the energy dependence is not completely flat, but it is smooth, and 
slow enough that its effect may be treated as a small correction. 
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dt 1 GeV 2 
[t = t min ] 
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FIG. 12: We obtain this fit using data from Ref. 



27 



3,3 



Ey (Gev; 



(listed in Ref. |l| ) . We use an inverse- 



fourth function at low energies and apply a least-squares fit. The peak in the parameterization 
yields a small effect from Fermi motion (at a few GeVs) because of the small mass of the p°. 

APPENDIX B: DEUTERON POLARIZATION 



In this appendix, we give an over view of the treatment of deuteron spin as it is presented 



m 



2l| . In order to evaluate the cross section, we must determine how the operator, ty* m , (k 



|)\E' m (k) acts on the spin-1 ground state of the deuteron. The non-relativistic deuteron wave 
function in momentum space is written in terms of S and D states via the formula, 



* m (k)= u{k) - 8~V 2 w(k)S 12 \q,m) 



(Bl) 



where, 



u(k) 



1 



rdrj Q (kr)u(r) 



V27T Jo 
1 f°° 

w(k) = / rdrj 2 (kr)w(r). 

V 2ir Jo 



(B2) 



The real functions, u(r) and w(r), are taken from any realistic model of the deuteron wave 
function, and in our computations we use the Paris potential |22j. The functions, j and j 2 , 
are the zeroth and second order spherical Bessel functions. In Eq. (jBip . \q, m) is a spin-one 
spinor representing the total angular momentum of the deuteron, and q is the quantization 
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FIG. 13: The low energy data here is from Ref. 16]. The curve at high energies was taken from 
Ref. |l|]. The dashed curve shows its extension to lower energies. The high energy data is taken 



from Ref. 



301 ] and is presented to establish the consistency of the high energy parameterization 



The curve at low energies has been fit to the low energy data of Ref. 



la] using a least-squares fit. 



axis. The tensor operator, Su acts upon the total angular momentum state to produce a 
sum over total spin states. In terms of the spins of the nucleons, it is given by: 

3 (a x ■ r) (a 2 ■ r) - ax • a 2 



S 



12 



(B3) 



The projection onto total spin states is 

S\0,q) 
S\-l,q) 

S\hq) 



yJ^-Y&OM - 1) - ^Y^m + ^Y 2 \6A)\l) 



16tt 



96tt 



y2{om-i) 



48tt, 



Y 2 -\9,<P)\0) 



96vr 



Y 2 ~\e^)\i) 



Yj{0A)\-i) 



^Y^m+J^Y^m- 



5 ^ ,T/I ' V 5 z v ,r/l ' V 5 

The functions, Y, are the usual spherical harmonic functions. With these equations, we 
can calculate the effective form factor for each polarization, and then sum and average over 
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final/initial deuteron polarizations. 
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